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Abstract 

We investigate the simulation of events with gaps between jets with a veto on 
additional radiation in the gap in Herwig++. We discover that the currently- 
used random treatment of radiation in the parton shower is generating some 
unphysical behaviour for wide-angle gluon emission in QCD 2 — )• 2 scatterings. 
We explore this behaviour quantitatively by making the same assumptions as 
the parton shower in the analytical calculation. We then modify the parton 
shower algorithm in order to correct the simulation of QCD radiation. 
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1 Introduction 



At modern colliders, such as the Large Hadron Collider and the Tevatron, a large fraction 
of the events produced are driven purely by the interactions of quantum chromodynamics 
(QCD). As a result of this, data generated in the early runs of the LHC is ideal for studying 
the nature of QCD. In addition to the availability of data, there have been many recent 
discoveries in the last ten years relating to the fundamental nature of QCD, such as 
non-global logarithms [Tl[2] and super-leading logarithms [3-[5]. From an experimental 
perspective, we also need a good understanding of QCD in order to reconstruct any 
physics of interest at said colliders. 

There are three fundamental theoretical problems associated with calculating observ- 
ables for purely QCD processes. They are confinement, the large coupling constant for 
small energies, and the large logarithms generated by radiative corrections. These three 
properties make fixed order calculations unreliable, since each order can be as important 
numerically as that which preceded it. Therefore the alternate method of resummation, 
which in some approximation contains all orders of effects, needs to be used. The resum- 
mation of large logarithms can be done in two ways, either analytically or numerically. 
The analytical approach calculates the first order corrections in the soft limit and expo- 
nentiates to obtain the leading term at each order of the all orders result. The numerical 
approach, referred to as a parton shower, evolves the process from the hard scattering 
down to some infra-red cut-off scale by randomly generating radiation according to a phys- 
ical principle. Formally the two approaches are equivalent in the large Nc limit, where 
Nc is the number of colours. There are a number of different parton shower programs 
currently available, see for example the recent review [6]. In this paper we will focus on 
the event generator Herwig-|--|- [IllH]. 

The gaps between jets process is an exclusive process which is ideal for working with 
early LHC data. By exclusive we mean that there is a restricted phase space for emissions 
of additional jets. In terms of the gaps between jets processes we impose a veto on 
radiation into the rapidity gap between the two hardest jets above a certain scale. From 
an experimental side there will be a large number of events produced since the underlying 
process is purely QCD 2 — )■ 2 scattering. The signal is fairly simple since all we are looking 
for is a number of jets of at least two. While we are vetoing additional jets inside the 
gap there is still a possibility of emission into the region outside of the gap. From the 
theoretical side the basic process is well understood at the lowest order of perturbation 
theory, often referred to as the Born level, and approachable at all orders using a leading 
logarithmic resummation. By varying the size of the rapidity gap and the amount of 
hadronic activity allowed it is possible to investigate a wide variety of QCD properties. 
In addition to this we can use this process to test how the parton shower, which is formally 
defined in the collinear limit, fares with regard to wide-angle radiation. Early LHC data 
on dijets with a veto have already been analysed by ATLAS with center of momentum 
energy y/s=7TeV [H]. 

The outline of the paper is as follows. First we summarise the current analytical 
approach to jet vetoing in QCD 2 — 2 scatterings in section two, as described in Ref. (TO]. 
Next, in section three, we explain how the analytical calculation differs from the parton 
shower picture in Herwig++. We then, in section four, explain how it is possible to 
modify the analytical calculation in order to quantitatively understand the parton shower 
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behaviour. Moving on to section five, we use the insight gained from the analytical 
approach to modify the parton shower so that the unphysical behaviour is removed. 
Finally, we state and interpret our results in section six and then give a conclusion 
on the status of the parton shower in section seven. While our focus is on Herwig++, 
we mention in several places the similarities and differences with respect to other parton 
shower algorithms. 



2 Gaps between jets — Analytical 

In this section, we summarise the work of Refs. p^HT3] . For the purpose of this paper we 
are mainly considering the colour structure, and as such will omit the details regarding 
the kinematics and convolution with PDFs. We will consider a partonic QCD 2 — )■ 2 
scattering. 



where each pi is an arbitrary parton. We require that there are two high transverse 
momentum hard jets separated by a rapidity gap Y with limited hadronic activity. An 
example t-channel process is shown in Figure [U We will refer to this as a gaps between 
jets event. 



Proton 1 



Proton 2 




Jet 1 



Rapidity gap 



Jet 2 



Figure 1: Dijet production with a rapidity gap. The black lines represent the partonic 
interaction and those with arrows can be quarks, antiquarks or gluons. The coloured 
particles are generated by radiative corrections. 



We restrict ourselves to the case of zero additional hard jets outside the gap. To 
quantify limited hadronic activity, we require that there are no jets in the gap above a 
scale Qo, which we refer to as the veto scale. We will also define an additional scale, Q, 
which is the average transverse momentum of the two leading jets. The act of limiting 
the amount of radiation inside the gap induces a miscancellation between the real and 
virtual corrections which generates a large logarithm of the form 

(2) 



Log 



Qo 
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In this paper we will consider Q to be some energy scale accessible by the LHC, and Qq to 
be small relative to Q but large enough that perturbation theory can be used. Note that 
since we are considering radiation from the hard jets into the gap, we will only generate 
single soft logarithms. To perform the resummation, it will be useful to write the matrix 
elements in terms of a basis of colour states. The partonic cross section at the Born level 
can be written in the form 

ao = Tr[HS], (3) 

where H is the hard matrix containing only Lorentz structure and S is the colour metric 
containing information about the interferences of the different colour flows. After the 
resummation of soft logarithms inside the gapEl we can write the partonic gap cross 
section as 

crGap = Tr[He-«r'Se-«r], (4) 

where F is the anomalous dimension which describes the evolution of the colour structure 
upon attaching an unresolved soft gluon and ^ is given by 

^ Jqo 

The exponential factors serve as a Sudakov form factor [13]. If we know H, S and F 
then we have determined the observable completely. The bulk of the work is that these 
objects are M-dimension matrices in colour space, where M can be up to six for general 
SU(A^c) theory, and if we wish to retain the full colour structure then each of the elements 
is a general function of A'^c- 

In order to calculate these three objects we need to choose a specific basis for the colour 
states. For the analytical calculation the choice of an orthonormal basis is convenient as 
the colour metric becomes the unit matrix. This was the basis used by Ref. [TD], and is 
similar to the basis used in Ref. pU] , but that basis was not normalized and did not benefit 
from the symmetric nature of the anomalous dimension in an orthonormal basis [15j. An 
alternate choice of basis, more suited to parton showers, which is the colour flow basis, 
will be discussed in the next section. 

The hard matrix is determined by the Lorentz algebra of partial amplitudes |16] . 
Specifically, if we write the matrix element Mj of the Feynman diagram i as an expansion 
on a basis of colour states | Cj), 

M, = Y^Af\C,), (6) 

3 

then the hard matrix is 

k,l 

■^Here we are only considering the resummation of the logarithms which occur due to the miscancella- 
tion between the real and virtual emissions inside the gap. There are additional non-global logarithms [T] 
resulting from secondary emissions into the gap which have not been resummed. 
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where the sums run over all possible Feynman diagrams contributing at the Born level. 
The soft metric can be determined by the projection, 

Sjj = (Cj I Cj). (8) 

By decomposing the colour basis states into Feynman diagrams, the components of the 
metric can be calculated by colour trace methods. The final task is to calculate the 
anomalous dimension, which can be written in colour basis independent notation as 

r = irt? + ^7rti-t2 + ip(t2 + t2), (9) 

where t, is a colour charge operator which quantifies the action of attaching a gluon to 
the parton i in colour space. The form of tj depends on the parton which is connected 
to the gluon and whether it is initial or final state. The first term, which is proportional 
to the size of the rapidity gap Y and the colour exchange across the gap t^, is due to the 
emission of wide-angle gluons. This term serves as the main suppression factor for the 
resummed cross section. The second term is due to Coulomb gluons which appear in the 
virtual corrections that connect partons that are either both in the initial state or both in 
the final state. This term is mainly responsible for mixing between different colour basis 
states. The last term is an edge effect where small-angle radiation at the edge of a jet 
is able to make it into the gap. This term provides an additional small suppression. For 
large values oiY (Y > 3) and for a jet radius oi R = 0.4 then p ~ 0.6 and independent 
of Y . In Figured] the wide-angle gluons are purple, the Coulomb gluons are blue and the 
gluons which generate the p factor are green. 

The calculation of the anomalous dimension can be done in a similar way to the colour 
metric. The anomalous dimension in an arbitrary colour basis is 

r., = E^rfc'(c.|r|c,). (10) 

k 

To compare to the experiment data, the resummed cross section needs to be convoluted 
with the PDFs. For a convenient choice of variables it is possible to factorise out the 
PDFs from the hard scattering and we are able to write the resummed hadronic cross 
section in the form, 

aH{Q,Q,,Y)=Y,^^^{Q.y>GapAQ^Qo,Y), (11) 

i 

where the summation is over all possible partonic processes and the function £j contains 
the PDF integrals. In the eikonal limit the emission of gluons does not change the 
kinematics and therefore the PDFs are independent of Qq even after the resummation. 

3 Gaps between jets — Parton Shower 

The Herwig+-|- parton shower, hereafter referred to as just the parton shower, is a se- 
quence of quasi-collinear 1 — )■ 2 splittings governed by the DGLAP equation [T7HT9] . The 
shower works by evolving a scale (f [20] , related to the angle of the splitting. The DGLAP 
equations are solved using a random number generator in order to generate the scale of 
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the next splitting, and the fraction of momenta to pass onto the two child partons. Here 
we will only consider the evolution scale and the colour structure of the parton shower. 

The colour structure of the parton shower differs from the analytical calculation in 
a number of ways. Before the parton shower begins, a hard process is generated in 
Herwig++ and a colour structure is chosen and retained throughout the whole evolution. 
The colour structures implemented in Herwig++ are a subset of the colour flow basis, 
where each basis state is a unique colour flow. In this basis, the colour metric is not 
orthonormal and thus there are interferences between colour states even at the Born 
level. To remove these interferences, the hard scattering is redefined in such a way that 
the Born level cross section is still generated [2T], but with the colour metric now being 
the unit matrix. Since the parton shower is probabilistic, the anomalous dimension is 
also taken to be real and diagonal. Thus we omit terms which are due to Coulomb 
gluons. This is equivalent to taking the large Nc limit in the subsequent evolutior0, since 
interferences generated by Coulomb gluons will in general be suppressed by factors of 
(1/iV^) relative to the leading contributions. Since the singlet colour structures have no 
associated initial hard scattering element, they will never contribute to physical processes 
in the large Nc limit. Hence these singlet colour structures are not implemented in the 
parton shower. Note that though the parton shower works mainly in the large Nc limit, 
there are also sub-leading corrections implemented in the colour charges of gluons and 
quarks, in order to generate the correct amount of radiation in the combined soft and 
collinear limit. 

Most of the preceding discussion applies to all other general purpose event generators. 
Hard processes are generated exactly and then projected probabilitistically onto the colour 
flow basis (or another that is equivalent to it in the large Nc limit) for subsequent evolution. 
However, in the next step, the evolution itself, one may draw a distinction between 
algorithms based on evolution of individual partons ("parton showers" like HERWIG, 
Herwig++ and the original virtuality-ordered shower in PythiaG [22]) or individual colour 
lines ("dipole showers", like the pt-ordered shower in PythiaG and PythiaS p^, Sherpa 
[2A\ . Matchbox [25], Vincia [26] and most other recent implementations). 

In the large Nc limit we can decompose the colour structure of a gluon into that of a 
quark line and an antiquark line. Any gluon produced in either the hard process or the 
subsequent shower will be connected by colour lines to two other partons. The scale of 
the radiation is proportional to the angle between the parton and the colour connected 
partner. For quarks there is only one possible partner, but for gluons there are two 
possible choices of partner. If we consider the case where one partner has a colour line 
scattered by almost tt and the other by a small angle, as occurs in, for example qg ^ qg, 
then there is a large difference in the two possible initial evolution scales. An illustration 
of the two different radiation patterns is shown in Figure [21 The current implementation 
in all three parton shower algorithms just mentioned is that the partner of a gluon is 
chosen from the two possibilities with 50-50 chance. This implementation generates the 
correct amount of radiation for inclusive events. The radiation in exclusive events is, 
however, not treated correctly. For jet vetoing processes the amount of radiation into the 

''Coulomb gluons are not the only source of sub-leading colour interferences. There are also interfer- 
ences between colour states suppressed by rapidity that exist even when Coulomb gluons are removed. 
However we have found that the difference between removing Coulomb gluons and taking the large Nc 
limit is minimal. 
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gap is highly dependent on which of the two partners is chosen. 




Figure 2: Left: Wide-angle radiation is produced with a colour factor Ca when the 
partner with the largest scale is chosen. Right: Only small- angle radiation is produced 
when the smaller scale is chosen. The blue and red lines are separated by wide and small 
angles respectively. The bold lines are the colour lines of the gluon and the faint lines 
are the colour lines of the partners. 



Figure 3: Colour line representation oi g gg splitting in the large Nc limit. 

Once the partner is chosen the initial scale of radiation is determined and the par- 
ton shower begins to generate radiation. Each time a splitting occurs the colour line 
connections must be generated. For q qg splitting there is only one possible colour 
connection. The parton shower picture for g ^ gg splitting is given in Figure El There 
are two possible colour line connections, which can be thought of as either the colour line 
or anti-colour line splitting individually. The parton shower chooses which line to split 
with a 50-50 choice. At the level of the shower this choice does not have any physical 
consequence. It is only when non-perturbative hadronization occurs that attaching ra- 
diation randomly can result in a different distribution of hadrons in the final state. By 
contrast, dipole showers evolve the two colour lines independently and the 50-50 choices 
of colour partner and emissions colour structure are not needed. 

In addition to the above changes there are also additional effects implemented in the 
parton shower that are not present in the analytical calculation. Some non-global loga- 
rithmic effects are included in the parton shower as a result of the angular ordering [27]. 
Energy-momentum conservation is absent from the analytical calculation due to the use 
of the eikonal approximation, but it is included in the parton shower. The implementa- 
tion of energy-momentum conservation into the analytical calculation is discussed at the 
end of the next section. 




Large 
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4 Modifications to the analytical calculation 



To understand quantitatively how the parton shower currently behaves we can construct 
an analytical model of the colour evolution. This amounts to modifying our current 
analytical calculation in a number of ways. First we define a new hard scattering matrix 
H such that its trace is equal to the Born cross section |2T] , in the following way 

Tr[HS] ^ Tr[H], (12) 



H,, = H,,Tr[HS]/Tr[H], (13) 

where we have normalized the colour basis states such that the diagonal elements of S 
are unity. We then take the large Nc limit in the anomalous dimension and drop any 
imaginary terms. Since the anomalous dimension is now diagonal it can be treated as a 
rescaling factor for the hard structure of each of the independent colour flows. We can 
write the diagonal anomalous dimension in the form@ 

r = + ^ + ^ ^ e"^ = e'^'e- ^« ^"e^ ^« , (14) 
g a 

where we have factored out the wide-angle radiation from the small-angle radiation, and 
further factored out the contributions to the wide-angle radiation from the different types 
of parton and the sums over i contain the contributions of all the partons of type i. 
For four external partons we will always have an even number of quarks and gluons 
separately and so we can treat the colour evolution of each parton type together. For 
the contributions to the anomalous dimension we only need to consider the number of 
colour lines crossing the gap. If a parton's colour line crosses the gap then it contributes 
a factor 



r. = ia, (15) 

where is Cp for a quark line and ^Ca for a gluon line. We implement the random 
choice of the partner in the analytical calculation by modifying the colour evolution factor. 
Since the parton shower generates colour evolution at the cross section level rather than 
at the amplitude level, and takes the anomalous dimensions to be diagonal, then we can 
simply double the anomalous dimension and take this to be our total colour evolution 
factor. At leading order in the anomalous dimension we can write the colour evolution 
factor for a single gluon as 

e-2r. ^l(l + e-^^0, (16) 

which is the mathematical equivalent to the statement that half the time we have no 
wide-angle radiation and half the time we have twice as much. In the case where both 
partners are on the same side of the gap then the left-hand side of equation f lTB]) is used, 
but when we have one partner on each side of the gap then the right-hand side of equation 
(|T6i) is used. 

^Here we only consider the colour structure and absorb the logarithmic factor into the definition of F. 
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From equation f[T6l) we can see that processes with hard gluons will tend to have 
a larger number of gap events than expected. For gg — )■ gg there is a possibility that 
the choice of colour will be a quasi- singlet, where all four choose partners on the same 
side of the gap. We refer to this as a quasi-singlet since it has non-zero Born level hard 
components, whereas a physical singlet can only gain a hard component by radiative 
corrections interfering two colour flows. The quasi-singlet arises from the treatment of 
radiation in the parton shower, not because of the group properties of the theory. For 
large values of L or F, the suppressive factor then becomes 



_ 2asCAP 

.--(HK'^'r-ie— i(|V - . (17) 



which scales like a power law, where the last approximate equality holds only for a fixed 
coupling Ug. In the analytical calculation the suppressive factor always has the form 



2asCA{p+AiY) 

Qo 



where Ai are constants greater than zero and i runs over all different colour flows. For 
large rapidity gaps this expression vanishes more quickly than that of the quasi-singlet. 
This quasi-singlet therefore results in a greater gap cross section at the more extreme 
edges of phase space. 

To illustrate the effects of these unphysical singlet type terms, we will look at the 
magnitudes of the different partonic processes in the resummation, with and without the 
modifications to the colour evolution. We define the resummed partonic fraction Ui of a 
process i to be, 

00. = (19) 

where ai^r is the resummed cross section for a given partonic process i. 

The upper plot in figure H] shows the partonic fractions for the analytical resummation 
without the changes to the colour evolution, for Y=5, (5=500 GeV and A/i=14TeV. We 
see that it is mainly quark-quark interactions that dominate. There are two reasons for 
this, the first of which is that we are in the high x regime (x ~ 1) where the quark 
PDFs dominate over the gluon PDFs. The second reason is that gluons have a greater 
colour charge and therefore tend to radiate more, which results in less events involving 
gluons that contain an empty gap region. The plot is cut off at Qo=l GeV as beyond 
that point it is mainly the non-perturbative behaviour which dominates and this is not 
implemented in the analytical calculation. The lower plot in Figure H] shows the partonic 
fractions for the analytical resummation with the changes to the colour evolution, again 
for Y=5, Q=500 GeV. We see that there is a large difference here, and that now the 
processes involving gluons begin to dominate for low veto scales. It is this quasi-singlet 
behaviour which is generating the large partonic fraction for the gg — ?■ gg process. 

In addition to the change in colour structure we also need to take into account ef- 
fects of energy- momentum conservation. In the eikonal approximation one assumes that 
gluons which are emitted from the hard partons are infinitely soft and do not change 
the kinematics of the Born process. This correctly produces the leading logarithms we 
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Figure 4: Partonic fractions of the resummed cross section as a function of veto scale for 
Y=5, (5=500 GeV. The upper and lower plots are without and with the change to colour 
evolution respectively. The PDF sets used are MSTWO8I0 



require for the resummation. In reality any emitted gluon carries a finite energy and mo- 
mentum and this should be taken into account. This affects the PDFs, since more energy 
is needed in the hard scattering, and the hard partonic cross section, since recoil adjusts 
the kinematics. To estimate the effects of energy-momentum conservation we consider 
the case of a single emission of a gluon of energy kt and rapidity y'. This modifies the 
arguments of the PDFs to be 
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Xl,2 = Xi,2 H 1= , (20) 



where xi^2 are the arguments of the PDFs without imposing energy-momentum conserva- 
tion, ^/s is the hadronic CMS energy and y is the sum of the rapidities of the two hard 
partons. We are no longer able to factorise out the and y' integrals previously contained 
in the anomalous dimension from the PDFs. Instead of this we consider the simpler case 
where kt=Qo and y'=Y/2, which will give an approximation of the maximum effect of a 
single emission. We then perform the calculation as before with the factorization, only 
now with the overestimate for the PDF arguments 

^i,2 = a;i,2H 1= — • (21) 



The effects of implementing partial energy-momentum conservation in the analytical cal- 
culation on the gg — )■ gg gap fraction are shown in Figure |3 We define the gap fraction to 
be the ratio of events with a gap to the total number of events. We see that after imple- 
menting energy-momentum conservation there is a much larger gap fraction for lower Qq. 
This is because higher energy emissions are suppressed by increasing x{Qo) in the PDFs. 
Note that while we have overestimated the effects due to one emission, we have neglected 
the effects due to multiple gluon emissions, and in some sense have underestimated the 
effects of energy-momentum conservation. 



5 Modifications to the parton shower 

Using the insight we have gained from the analytical calculation we now understand how 
the parton shower behaves relative to what we would have expected. The current parton 
shower generates an excess of wide-angle radiation 50% of the time and no wide-angle 
radiation the other 50% of the time. This leads to a larger number of events passing the 
veto on radiation than one would expect. To correct this behaviour we choose to modify 
the internal structure of the parton shower. 

Our aim in changing the parton shower is to allow each line of a gluon to radiate 
quasi-independently i.e. in a more similar way to the dipole shower formulation. The 
first modification is to change the assignment of partners after the hard scattering so 
that the partner with the largest angle is always chosen. When this is done the shower 
records the scale with respect to the furthest parton, qj, and the scale with respect to 
the other parton, g„. The shower then begins with half the colour factor, ^Ca, as only 
one of the colour lines can radiate at the largest angle. The colour structure is set up 
such that the emitted gluons are only attached to the wide-angle colour line. Once the 
scale has been evolved down to g < g„ then both lines are able to radiate and the colour 
factor is restored to Ca. The shower will attach the colour lines of the additional gluon 
with a 50-50 chance to one of the two lines of the parent. 

The radiation pattern for the modified parton shower is shown in Figure El Note that 
the emitted gluons are connected to the correct lines, which results in a change to the 
hadronization behaviour. 
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gg->gg Gap Fractions (Q = 500 GeV, Y = 3, S = (14 TeV) ^ ) 




QO / GeV 

Figure 5: Effects of partial energy-momentum conservation in tlie analytical calculation 
on the gg — > gg gap fraction. The implementation of energy-conservation increases the 
gap fraction for low Qo? as observed in the transition from black to red and green to blue 
respectively. 




Figure 6: The radiation pattern for the modified Herwig-|-+ shower. The blue and red 
lines are separated by wide and small angles respectively. The bold lines are the colour 
lines of the gluon and the faint lines are the colour lines of the partners. 

6 Results 

For our investigation we choose to use the MSTWO8I0 PDF set [21] and use SISCONE p] 
via Fast Jet |30j with cone radius i?=0.4 and overlap parameter 0.5 for jet finding, as was 
done in the previous analysis. For the first two sets of results we turn off hadronization 
and multiple interactions, since these are not simulated in the analytical approach. 
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Now that the corrections to the parton shower have been made we will validate the 
approach by considering the gg — )■ gg gap fraction, which is the most sensitive to the way 
radiation is treated in the parton shower. In the absense of energy-momentum corrections 
the gap fraction is simply the ratio of the resummed cross section to the Born cross section. 
There are three choices of variable that we can use in order to investigate the gap fraction, 
the hard scale Q, the rapidity gap size Y and the veto scale Qq. Changing Q or Y will 
modify the PDFs, the hard components and the colour suppression due to the large 
logarithms. In order to investigate the effects due to the change in the colour structure 
only, it is therefore better to investigate the gap fraction using Qq, which changes only 
the large logarithms and, more weakly, the PDFs. Using Qq as a variable has previously 
been suggested for other processes 



gg->gg Gap Fractions (Q = 500 GeV, Y = 3, S = (14 TeV) ^ ) 
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Herwig++ 2.5.0 with the colour evolution modification 
Analytical model of the parton shower 

Analytical model of the parton shower with the colour evolution modification 
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Figure 7: The gap fraction as a function of Qq. The red and black curves are the 
numerical results with and without the modifications to the colour evolution. The green 
and blue curves are the analytical results with and without the modifications to the colour 
evolution, both including the modifications to account for energy-momentum conservation 
described in Section 4. 



Having decided upon the choice of variable, we now plot the gap fraction as a function 
of Qq in Figure [71 We choose the parameters 1^=3, (5=500 GeV as this enhances the 
logarithmic wide-angle suppressions and take the beam energy to be y^=14TeV. It is 
clear that none of the curves seem to match and the analytical calculation generates more 
high energy radiation than Herwig++. Since we are interested in validating our change 
to the parton shower and not the gap fraction itself, we can remove any additional 
differences between the numerical and analytical approaches by forming ratios. Two 
different methods of dividing out are shown in Figures |8] and El From the ratio plots we 
can conclude that, although the analytical curves do not match the numerical curves, our 
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modification performs as expected. Tliere are some minor disagreements in tlie low Qq 
region but since tliis region is dominated by non-perturbative effects we can no longer 
trust tlie analytical calculation, and the behaviour in the parton shower begins to be 
influenced by the infra-red cut-off. 

gg->gg Gap Fraction Ratios (Q = 500 GeV, Y = 3, S = (14 TeV) ^ ) 




Herwig++ 2.5.0 with the colour evolution modification divided 
by Herwig++ 2.5.0 

Analytical model of the parton shower with the colour evolution 
modification divided by the analytical model without the modification 
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Figure 8: The ratio of the colour evolution modified gg — )■ gg gap fractions to those 
without modifications for Herwig++ 2.5.0 (black) and the analytical model of the parton 
shower (red). 



Having convinced ourselves that the modifications have the correct physical behaviour, 
we now move onto the gap cross section as a function of Q, which was the original 
motivation for our present analysis. This is shown in Figure [TOl The upper plot shows 
the resummed gap cross section, and the lower plot shows the ratio of the gap cross 
sections to the full colour structure result of Forshaw et al. [10]. The analytical predictions 
for the parton shower with and without our modification are clearly below that of the 
FKM result throughout the whole range of Q. Again we see a difference between the 
numerical and analytical results, with the numerical results always being greater than 
their analytical counterparts. Performing the same ratio analysis as above we have found 
that the modification behaves as expected. 

So far we have only considered the change in behaviour after running the parton 
shower. Since partons are not observed in the final state we should also include the effects 
of hadronization. These are shown in Figure [Til We see that including the hadronization 
effects change the gap fraction by a small amount relative to the parton shower level, with 
noticeable effects in the low Qq region. We do expect a difference because, comparing 
Figures [2] and [HI we see that the current algorithm has the effect of "pulling" the colour 
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gg->gg Gap Fraction Ratios (Q = 500 GeV, Y = 3, S = (14 TeV) ^ ) 
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Herwig++ 2.5.0 divided by tlie analytical model of the parton shower 

Herwig++ 2.5.0 with the colour evolution modification divided by 
the analytical model of the parton shower with the same modification 
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Figure 9: The ratio of Herwig++ 2.5.0 to the analytical model of the parton shower with 
the colour evolution modification (black) and without (red). 



line from the small-angle scattering out into the gap region, which does not occur in the 
modified algorithm that we have presented. 



7 Conclusion 

Herwig++ is very successful at making a correspondence between hard processes and 
experimental observables. This is mostly due to the parton shower, which handles the 
evolution of a hard 2 — )■ 2 scattering into experimentally observed high multiplicity final 
states. The underlying mechanisms of the parton shower are obtained from our under- 
standing of the theoretical nature of QCD. As such, it is important that we investigate 
and correct any possible deficiencies in the parton shower approach. 

In this paper we have shown that the current parton shower of Herwig++ did not 
produce the right amount of wide-angle radiation for inclusive events. We interpreted 
this behaviour using an analytical model of the parton shower behaviour in gaps between 
jets events. This led to a modification of the internal structure of the parton shower. 
Although the changes we have made may have small effects in the current phase space 
the physical behaviour of the parton shower is now more in line with the analytical picture 
of radiation arising from QCD. 

We would expect similar conclusions to apply to other parton shower algorithms 
(HERWIG and Pythia6), but more recent algorithms based on a dipole picture already 
include the effect we have considered. 
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Figure 10: The gap cross section for all possible QCD 2—7-2 processes as a function of 
the hard scale Q. 
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Figure 11: The ratio of the gg — )■ gg gap fraction generated by Herwig++ with hadroniza- 
tion to that generated by Herwig++ without hadronization. 
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Whether the parton shower is still successful with regards to predictions in gaps 
between jets events at the LHC has yet to be determined. Both the numerical and 
analytical approaches have deficencies with regards to each other. The current level of 
data analyses at the LHC [9J is not enough to distinguish between the two approaches. 
As more data arrive we will be able to see which of the two approaches best describes 
nature. 

From a theoretical point of view however, the parton shower is still lacking the inclu- 
sion of physical singlets responsible for a significant fraction of the cross section at large 
Y and Q. These are not included since they involve evolution at the amplitude level, 
which is not possible as the current parton shower evolves the hard scattering at the 
cross section level. A full treatment of QCD in the parton shower will therefore require 
an amplitude level parton shower. 

The code to implement our modificatins has been released in Herwig++ version 2.5.2, 
and may be switched on using the instructions in the Appendix. 
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A Switches 

To enable the new colour evolution model in Herwig++ 2.5.2, add the following lines to 
the "..in" file used with Herwig++ 

cd /Herwig/Shower 

set Evolver : ColourEvolutionMethod 1 

set PartnerFinder : PartnerMethod 1 

set GtoGGSplitFn:SplittingColourMethod 1 
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